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Abstract 

As generalizations of the fermion seniority model, four multi-mode Hamiltonians are 
considered to investigate some of the consequences of the pairing of parafermions of order 
two. 2-particle and 4-particle states are explicitly constructed for = —G A^A with 
A"!" = ^J2m>o ^In^^-m ^"^^ t^c distinct He = —G G'^C with C"^ = \ Yl,m>a c^mC-mi and for 
the time-reversal invariant = -G (^t - C'^){A - G) and = -G {A^ + G'^){A + C), 
which has no analogue in the fermion case. The spectra and degeneracies are compared with 
those of the usual fermion seniority model. 
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1 Introduction 

The physics of fermion pairing and fermion condensates [1] is important in contemporary ele- 
mentary particle physics in precision QCD calculations for hadron spectroscopy (e.g. via lattice 
gauge theory or chiral effective Lagrangians) and in research on dynamical electro- weak symmetry 
breaking of the standard model (e.g. via technicolor or a tt condensate). In this paper, we con- 
struct and study four Hamiltonians as generalizations of the fermion seniority model [2] in order 
to investigate some of the consequences associated with the pairing of parafermions of order 2. A 
physical significance of "order p = 2" is that 2 or less such parafermions can occur in the same 
quantum state. Usual fermions correspond to p = 1. 

Although the idea of the possible existence of fundamental particles associated with other rep- 
resentations of the permutation group is an old and simple one [3] , and despite the existence of 
significant general results in relativistic local quantum field theory concerning properties of ele- 
mentary particles obeying parastatistics [4,5,6], calculations in this field are sometimes intractable 
because of algebraic complexities arising from the basic tri-linear commutation relations, see (1-3), 
versus the standard bi- linear commutation relations which occur in order p — 1. Order p = 2 is 
indeed simpler than p > 2 because there is a "self-contained set" of 3 relations [4], see (7-9). 
With respect to representations of the permutation group, consideration of the "order p — 2" 
parafermions is not special in that there are still the two d = 2 dimensional representations, 
a mixed representation still occurs ai d — 3, and multiple mixed representations still occur at 
d — A. On the other hand, for p > 2, mixed representation(s) with both totally-symmetric and 
totally-anti-symmetric ones do occur, starting at d = 3. 

The first Hamiltonian considered in this paper is Ha = —G A^A with A = J2m>o where 
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j^{m) _ lc_^Cm- The mode index k, 1,171 ranges from 1 to O and J2m>o denotes summation over 
1,2, ...Q. In this paper, summation symbols are always displayed, so repeated indices are not to 
be understood to be summed. In many treatments of the usual p — 1 seniority model, for instance 
in nuclear physics applications, fl is the number of {I,— I) pairs and the "mode index" k,l,m is 
the magnetic quantum number. In this paper, the time-reversal operation, T, will be analogously 
defined to exchange I —I , but except for the use of this exchange operation, no explicit 
physical significance such as "magnetic quantum number" is associated with the k,l,m index. 

In paraquantization, it is instructive to begin by summarizing the parafermi and parabose 
cases together: The basic commuation relations are 

[Cfe, [c},cj=p] = 2 Ski Cm, (1) 

[ck,[clcU^] = 2SkiclT'2Skm4, (2) 
[cfc, [Q,c„]zp] = (3) 

following the standard convention that the upper (lower) signs correspond respectively to 
parafermions (parabosons) . The minus subscript is often suppressed so [A, B] = [A, -S]- = AB — 
BA, and {A, B} = [A, -B]+ = AB+BA. In this paper, the corresponding creation and annihilation 
operators for the ordinary p = 1 fermions are labeled al and a^. For p = I, a]^aL^ = — aLjn(4n ' ^ut 
for p > 1, cl^c^m ^-nd c^-mC^m distinct operators. The number operator for the para-particles 
is Nk = |[Cfe,Cfc]:p ± I with the order of the paraparticles, p, being the maximum number of 
parafermions (parabosons) in a totally symmetric state (anti-symmetric state). We assume a 
unique vacuum state with the usual properties Cfe|0 >= 0, < 0|0 >= 1, and CfeC; |0 >= p Ski\^ >■ 
From here on in this paper, the cj, and Cm are parafermi operators of order 2. The following 
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two commuting operators frequently occur in this pairing analysis: 



N=-Y1 ([4, Cm] + [elm, C-m]) + 2 Q 



(4) 



m>0 



P=nYl ({4> Cm} - {elm, C-m}) 



(5) 



m>0 



Note that N is the sum of the parafermion number operators for the 20 modes. However, although 
c|, and Cm are parafermi operators of order 2, p has the formal structure of being the difference 
of parabosonic number operators for the m > and m < modes. The appearance of this T- 
odd operator p is a noteworthy difference versus the ordinary p = 1 seniority model, in which it 
vanishes. 

Since {B^"^^y — |c|„cL^ and {D^"^^y — ^clj^cl^ are distinct operators, we also consider a 
second Hamiltonian He = —G C^^C with C = J2m>o D^"^\ In A and in C the parafermions with 
m and —m are paired and so states constructed as polynomials of A and C will be labeled as 
states of seniority zero, s = 0, since they are states built out of paired particles. Thus, for Ha,c, 
the seniority s is the number of unpaired parafermions in the state, just as in the p — 1 case. 

In the investigation and analysis of fermion pairing phenomena in condensed matter physics, 
the presence or absence of time reversal invariance and its consequences has been one of the im- 
portant symmetry considerations [7]. If, in this respect, one does treat the m index on Cm as 
corresponding to a magnetic quantum number, then under the time reversal operation, B^"^^ ^ 
D^"*) and Ha ^ He ■ Although there appears to be no obvious violation of time-reversal in- 
variance at the observable's level corresponding to such a discrete switch of Hamiltonians, e.g. 
spectra will be the same for the two respective Hamiltonians, we find this a somewhat radical 
formal situation which appears not very easily generalized, in particular with respect to inclusion 
of kinetic energy terms and perturbations. Accordingly, in this paper we also consider two time- 
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reversal-invariant Hamiltonians which are mapped into themselves by this time-reversal operation. 
These are = -G {A^ - Ct)(A - C) and H^+) = -G {A^ + C^){A + C). Note that does 
not exist in the p — 1 case. As will be discussed below, results for such as its spectrum can 
normally be obtained from those for by appropriate "substitution rules." 

The operators N and p commute with Ha,c but only N, not p, commutes with H(^:^y 
As in the usual quasi-spin formalism [8] , Section 2 of this paper treats the algebras associated 
with the Ha,c Hamiltonians, N, p ,A, C, and other such two-body operators . Analogous to the 
s = operators A and C, two sets of s = 2 two-body operators Bi and Di with i— 1, Q — 1 are 
introduced. In Sections 3 and 4, these two-body operators A, C, Bi and are used to explicitly 
construct A/^-particle states with various senorities (N is even). These results are used to study 
the spectrum for H^ c and for if^p in comparison with that of the fermion seniority model H = 
~G '}Zm>o'An'^^-mHi>o'^-i'^i which has the A^— state spectrum 

E^{N) = -h3{N - s){2VLp=^- N - s + 2)-, s = 0,2,...A^. (6) 

The 2-particle and 4-particle states are exphcitly constructed for Ha,c in Section 3 and for 
in Section 4. In both cases for s = 4 the construction of the 4-particle states is only for 
Q: — A. For Ha,c and for ffqi, it is found that the necessary mutual orthogonahty properties of 
the 4-particle states are somewhat awkward to arrange using the two-body operators A, C, Bi 
and Di. There are two built-in "parafermi p-saturation" conditions for p = 2: (cj[.)^ = and 
(A^)^ = (C^)^ = when M = Q 1 [this second fact is also true in the p = 1 case]. For 
some of the additional 4-particle states, for instance, this has the consequence that some state 
normalization constants vanish for small Q values, because the states do not then exist. Some 
results for arbitrary VL beyond A = 4 are derived. In Section 3, in all cases, the spectrum and 
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degeneracies for Ha,c is found to be that of the usual p — 1 senority model. 

In Section 4, for for A'" = 2,4 by explicit construction of orthonormal states, a sizable 

number of additional states not present in the analysis of Ha are found to occur. For -f^(-), results 
are obtained for arbitrary N-particles states which can be constructed as polynomials in only the 

and operators. In all cases, the spectrum of is found to be that of the p—1 senority 
model, except that flp=i in (6) is replaced by 2fl. However, for there are many additional 
degeneracies beyond those of the usual p—1 senority model. These degeneracies can be specified 
by an appropriate use of the seniority number, s. 

A primary motivation for studying the seniority model is because it is a simple model which 
has been used for fermions to theoretically investigate and exhibit consequences of fermion-pairing, 
of the microscopic realization of superconductivity, and thereby of spontaneous symmetry break- 
ing. The most surprising result of this paper's analysis is that in this multi-mode framework for 
parafermions of order two, it is indeed possible to algebraically investigate the spectrum for each of 
the four Hamiltonians. In hindsight, this tractability is partially a consequence of three facts: (i) 
in a single mode, the pairing-operators B^'^^ = ^c-mCm and = ^CmC-m separately lead to a 
two-body operator, quasi-spin Hamiltonian structure which is similar to that of the fermionic case, 
(ii) for different modes, six Hermitian pairing-operators mutually commute: These operators are 
Bi^'") = + 52^"*) = - and Sg^"*) = and analogously 

for and eJ^^] = dJ^")] = dJ™^] = for / 7^ m where a,b=l, 2, 3, and (iii) 

for these operators in the same mode also [B^\ D^"^'^] = 0. 
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2 Two-body operator algebras 

For parafermions of order 2, one has the relations 

c\ciCm + CmClc\ = 2 5kl Cm (7) 
Ckclcjn + Cmc\ck = 2 5^ + 2 Cfe (8) 
CfeQCm + C^QCfe = (9) 

plus the Hermitian conjugate relations. Note from the left-hand side of these relations that there is 
a simple left -^-^ right reordering symmetry. On the vacuum state c^cj |0 >= 25fc;|0 >. A consequence 
of (9) is the "parafermi p-saturation" conditions noted in the "Introduction" . Useful commutators 
involving c|. , cl^, ... pairs are in the appendix of this paper. 

It follows that the s = two-body operator A = J2m>o B^^^ —\ Y^rn>Q C-mCm has the quasi-spin 
algebraic relations: 

[A, At] = -2Zas, [Za3, = A\ [Za3, A] = -A (10) 

where A^ = Ai + iA-2 with i = A3 = Zaz- For ^ = {Aif + (^2)^ + {ZasY, one finds 

[Sa ,^1,2,3] = and Ha = -G{Sa - {Za^Y + ^as)- Since [A^A^Za-a] = 0, Ha and Zas can 
be simultaneously diagonalized. The Ha eigenvalues are Ea — —G{sa{sa + 1) — {^sa )^ + z^a}- 
The explicit N-particle parafermion states of various seniorities corresponding to this spectrum 
are constructed in Section 3 below. In terms of the operators N and p defined by (4,5) 

Za3 = ^{N-p-2Q) (11) 

Similarly, for C = Em>o D^""^ = ^ J2m>o CmC-m, 

[C,C^]^-2Yc3, [Yc3,C^] = C^ = C, + iC2, [Yc3,C] = -C (12) 
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where 



yc3 = C3 = -(iv + p-2fi), 



(13) 



Since [C^C, 1^3] = 0, He and Ycz can be simultaneously diagonalized. The Hq eigenvalues are 
Ea = -G{sc{sc + 1) - ivzc Y + Vzc}- 



and 7V|0 >= p|0 >= 0. Unhke in the p — 1 case, here for the s — Q N-particle states due to the 



occurrence of p as well as iV, while N{A^)^\Q >= 2M(y4t)^|0 > and N{C^)^\0 >= 2M{C^)^\0 > 
for M a non-negative integer, there is a minus sign in p{C^)^\0 >— —2M{C^)^\0 > versus 
p{A^)^\0 >= 2M{A^)^\0 >. These states have respectively the energy eigenvalues E^'^''^\2M) = 
-GM{n- M + 1). 

There are the useful commutators 



Useful eigenvalues for the vacuum state are 




(14) 



[AT, At] 




(15) 



[Ha,A^ 



(16) 



(17) 



The p again appears with an extra minus sign in the analogous expressions 



[TV, Ct] 




(18) 



[He, Ct] 



(19) 



(20) 



Note [Ha,C^] = [Hc,A^]^0. 
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To explicitly construct states with seniority s 7^ 0, we define the s — 2 two-body operators 



B] ^ (21) 
= ^5Mt_^5(m)t. i = (22) 



m=l 



and 

D\ = ^(E ^^-A - ^ = 1' - 1) (23) 

^ i=i 

The 2-particle states S||0 > and dJ|0 > with s = 2 respectively have zero Ha^c energy eigenvalues, 

A'' eigenvalues of 2, and respectively p eigenvalues of ^2 like respectively A^C'^ . For VL arbitrary, 

[A, B\] = -2Z^Bi. [A\ Bl] = [El S]] = (24) 

[Ha^ bI] = 2GA^Z^Bi. [TV, Bl] = 2St, [p, sj] = -25] (25) 



where 



1 * 

^3Bi = 7(E{49 - C-jcLj} - i{c{.+i)C(i+l) - C_(i+i)cL(.+^}) (26) 



and 



[Z^Bi,A^] = 5l, [Z3^,5l] = Sl (27) 
[^3A, -^3Bi] = [-^asi, -^asj] — (28) 

On the vacuum state, Z^Bi\^ >— 0, so A5||0 >= Sj^l^jO >= 0. Alternatively, in terms of mode 
operators 

Z.Bi ^ \{j:{Nt~'-p^B^}-^{W'^-p%^'^}) (29) 

^ m=l 



= E^s"^-^^'^'^ (30) 

m=l 



where 



(m) 



B — 



Pb 



7M — 
^3 = 



^ ([cL> Cm] + [dm, C-m]) + 2 fl 



For f2 > 2, these operators do not completely close at the level of the A's,C"s,i?i'f 



instead involve mode operators: 



[Z3B^,B]] = { 



where i < denotes the smaller oi 
Similarly, 

[C, Dj] = -2Y^ni, [C\ Dj] = [D}, D]] = 
[He, Dj] = 2GC''Y,Di, [N, Dj] = 2DI [p, D}] = 2DI 



where 



1 ' 

Y3Di = j{^{cljC-j - Cjc]} - i{cL(i+i)C_(i+i) - C(i+i)c|.+i)}) 



4%i 



m=l 



and 



[>^3C,1^3Di] = [>3Di,l3D,] = 



On the vacuum state Y^Dil^ 0' so C£>J|0 >= AC^|0 >= 0, and 

YsDi ^ l{t{N^^^ + pP}-^{Nt'^ + f^^'^}) (42) 

^ m=l 

= ±Yt^-zYt'^ (43) 

m=l 

where Y^"^^ = | (^cLmC-m — Cmcln) ■ The states 5||0 >, A^|0 > have Zas eigenvalues of (1 — |n), 
and Yc3 eigenvalues of (— |f2); similarly dI\0 >, C^O > have Yes eigenvalues of (1 — ^Q), and Zas 
eigenvalues of (— |J1). 
Note that 

[Bi, B]] = -2Z3Bi< - + l){c{.^i)C(i+i) - c_(,+i)C^_(.^^)} (44) 

where again i < denotes the smaller of the last term's factor also appears in (26). On the 
vacuum state -Bj-S] |0 >= i{i + l)6ij\\0 >, so Bj\0 > are orthogonal for different i values. Similarly, 

[A, ^J] = -'2Y3Di< - ^Siji{i + l){cL(.^i)C_(i+i) - C(i+i)c{.+i)} (45) 

and A-DJ|0 >= i{i + l)Sij\\0 >. 



3 Spectrum of Ha 

Useful relations for treating arbitrary N particle states include: for r = 1, 2, ... 

A{A^y\Q > = r{n - r + l){A^y-'\Q > (46) 
{Ay'{A^Y\0 > =r!Q(Q-l)---(Q-r + l)|0> (47) 
ZA3(A^y\0 > ^(r-^n)(A^y\0> (48) 

Since we find the Ha case to be relatively simple, for instance it has the same spectrum as in the 
usual p — 1 senority model, we do not evaluate normalization constants in this section. 
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For arbitrary N >2, the s — states \N^ >o— {A^)^/'^\0 > have energy eigenvalues '{N) — 
-lGN{2n- N + 2). 

For the s = 2, N-particle states with N > 2, \N^ >2= st(At)(^-2)/2|o > have 4^^(Ar) = 
-\G{N -2){2Q- N). 

For s = 4, we are interested in both the Ha states for themselves and also for comparison 
below with those occurring in the analysis of the states: For fl — 4, the explicit orthogonal 
4-particle states include two with zero Ha eigenvalues. From the analogous p = 1 spectra and using 
completeness, we classify them as s = 4 states. These two states are |4^ >4— — S|)|0 > 

, |4^ >4= -^bIIbI + A^)\0 >; to achieve orthogonality these states are somewhat complicated 
in terms of the B} operators. The other Q = 4 states are special cases of ones discussed above: 
|4^ >o= (At)2|0 > , |4f >2= A^Bj\0 > where i = 1,2,3. 

Next, for the s = 4, N-particle states with iV > 4, we consider 
\N^ >4= S|5]At)(^-^)/2|0 > with i < j. Note that it is at the s = 4 senority level that in using 
the B} operators, the orthogonality requirements started to produce complications for the Ha case, 
and similarly in the case below. So in considering only \N^ >4, we are ignoring this difficulty. 
When Ha operates on this state, one can commute A past the first Bj using [A, B}] — —2Zj,Bi- 
In commuting Z-^Bi past the next ij] one can use (34), [Z^BijBj] = B\ since i < j, however, 
B} produces an (TV — 2)-particle state. We proceed by dropping such terms because such {N — 2) 
terms will not contribute if the Ha expectation value is calculated. In this manner, we obtain 
Ei'^\N) = -|G(A^ - 4)(2n - N -2). This is not a complete derivation because the \N^ >4 are 
not mutually orthogonal. If we proceed similarly, for arbitrary senority s — 2t, then for N-particle 
states with N > s, |iV^ >s= {-B^-B^^ • • • -B|j(At)(Af-^)/2|o > ^jth ii < ^2 < • • • < ^^, we obtain 
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E(^\N) = -\G{N - s){2n- N - s + 2); s = 0, 2, ...iV by induction. 



4 Analysis of 

For the Hamiltonians if(=p) whereas [A^, H(^zf)] = , 

[p, = tAG{A^C - C^A) = tA{Hac - Hca) (49) 

so the N-particle eigenstates of can no longer be classified by the eigenvalues of p. Note that 
H^^) = Ha± Hac ± Hca + He where Hac = G^tCt and Hca = GC^^t. 
The algebra associated with includes the equations 

[A±C,A^±C^ = -7Vac3, [^ + C,^t_ct]=p (50) 
[N,A^^C^ = 2(>lt^Ct)^ [p,At^Ct] = -2(^t±c't) (51) 

where Nac3 ^ (N - 2Q) = 2Za3 + '^Ycz = | Em>o ([c^j Cm] + [ct-^, c_^]) includes the zero point 
energy. In comparison, note that p = —2Za3 + 21^3. Also 

[F(_),At-Ct] = ^(At _ c^t)jv^^3 (52) 
[if(_),At + Ct] = -^(At-Ct)^ (53) 

A'^ = 2, 4 particle states: 

We list the orthogonal 2-particle states, their energy eigenvalues, and normalization constants 
N,{N) = s<N\N 
Ones with s = 0, 

|2- >o= (A'f - Ct)|0 >, E^~\2) = -2Gn] Noi2) = 2Q 
|2- >o= (A^ + Ct)|0 >, 4-^(2) = 0; No{2) = 20 
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Ones with s — 2, 

\2. >2= {B\ - Dl)|0 >, ^^^(2) = 0; N^{2) = 2^(^ + 1) 
|2r >2= {B\ + D\)\Q >, 4-)(2) = 0; N^{2) = 2i(i + 1) 
In calculating degeneracies, it is to be understood that i,j = 1, 2, ...,fl — 1. 
The additional states not present in the Ha,c cases are denoted with a "breve", or "short 
vowel" , accent. Such states are therefore not analogous to ones in the usual p — 1 seniority model. 
We similarly list the orthogonal 4-particle states: 
Ones with s — 0, 

|4- >o= (^t _ ct)2|o >^ E^Q-\4) = -2G{2Q - 1); No{4) = 2Q{4Q - 1) 
|4- >o= (At - Ct)(At + Ct)|0 >, Ei-\^) = -2^(0 - 1); 7Vo(4) = 2{^f 
|4r >o= ({At}2 + {Ct}2 + 2{^)A^C^)\0 >, Et\4^) = 0; iVo(4j = 4(20 - 1)(0 - 1) 
For instance in |47 >o, a "number subscript" is used on the additional states label to denote 
ones constructed with an ft dependence so as to achieve orthogonality. Due to the parafermi 
saturation such a state is absent for fl sufficiently small; this is seen in the norm vanishing 
and in the vanishing of a factor hke (^^)- Completeness in each "N,s sector" is achieved by 
construction. 

Ones with s = 2, which are orthogonal for i ^ j , 

|4- >2= {B} - Dj)(At - Ct)|0 >, 4-^(4) = -2G{2n - 2); Ar2(4) = 4(Q - + 1) 
|4r >2= {Bj + Dl){A^ - Ct)|0 >, ^-^l) = -2G(2f] - 2); N^i^) = 4(f^ - l)^(^ + 1) 

>2= {n{B}A^ + DjC^} + 2){i?lCt + dIA^)\0 >, 4-^(40 = 0; N^iA^,) = - 

i)(n-2)i(i + i) 

>2= (n{BlA^ - D\C^ + (O - 2){BlC^ - D\A^})\Q >, 4^^(42) = 0; NSi2) = ' 
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l)(0-2)i(i + l) 

Ones with s = 4, for the case = 4; all with zero energy eigenvalues and the same normaliza- 
tion constant A^4(4) = 144, 

|4- >4= {Bl - -Bl- {Dl - D|})|0 >; 

|4- >4= {b\ - d\){bI - b| + {d\ - d\})\q >; 

|4,- >4= (Sl + D\){Bl -Bl- {Dl - D\})\Q >; 

|4- >4= + d\){bI -bI + {Dl - Dim >; 

Plus two analogues of Ha states, and two additional states, 
\K >4= {Bl{Bl + At} + Dl{Dl + Ct})|0 >; Ar4(4j = 336 
|4- >4= + At} - Dl{Dl + Ct})|0 >; 7V4(4J = 336 

>4= 73(^24 + dIbI)\Q >; 7V4(4^) = 144 
|4^ >4= ^(5^4 - dIbI)\Q >; iV4(4^) = 144 

Some A/^-particle states constructed as polynomials in and Ct : 

Useful relations for treating arbitrary N particle states include: for r = 1, 2, ... 

(A-C)(At -Ct)'-|0 > =r(2Q-r + l)(At-Ct)'-i|0> (54) 

(A ^ C)'-(At ^ Ctj'-lO >= r! ^ r p(0 - 1) • • • (n - r + i + 1) 

f](f]-l)---(O-t+l)|0> 

and 

]V(At)'-i(Ct)'-2|0 >= 2(ri + r2)(At)^i(Ct)'-2|0 > (55) 
p(At)'-i(Ct)^2|0 >= -2(ri - r2)(At)'-i(Ct)^2|0 > (56) 
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Ignoring the mutual orthogonality requirement, for arbitrary N, the following seniority Sp=i — 
states are found to have the following associated eigenvalues: 

For > 2, {A^ - Ct)^/2|0 > has E^~\n) = -\GN{An -N + 2). 

For N>2, {A^ - Ct)(^-2)/2(^t + c^^\q > has E^o~\n) = -lG{N - 2){in - N). 

For TV' > 4, (At - Ct)(At + Ct)(^-2)/2|0 > also has Et\N') = -\G{N - 2)(4Q - N). So for 
N > 4, both these states are degenerate with the s — 2 states with N-particles. 

For n" > 4, (At - Ct)(^-4)/2({At}2 + {Ct}2 + 2{^}A^C^)\0 > has Et\N") = -\G{N - 
4)(4Q — N — 2) so this state is degenerate with the s — A state with N-particles. 

For these states, ignoring orthogonality, the spectrum is found to be the same as that of 
Ha, except that Jlp=i in (6) is replaced by 2Q. There are, however, many additional degeneracies 
which can be counted via completeness in each the N, Sp=i sector, for fl sufficiently large so that 
the absence of states due to p-saturation can be ignored. 

The states of the spectrum follow isomorphically by letting C — > — C and — > —D}: for 
instance, |2+ >o= (At + C;t)|o > has -^^+^(2) = -2Gf7; |2+ >o= (At - Ct)|0 > has E^-^\2) = 0; 
|2+ >2- {B} + dI)\0 > has 4^'(2) = ; and |2+ >2- {B} - dI)\0 > has 4^^(2) = 0. 
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Appendix: c^t, c_fct,... pair commutators 

Although equations (7-9) are sufficient, the following commutators are algebraically sometimes 



15 



more direct or useful as checks: 

[cfcC_fc,c|cLj = 2 5-k,i{^kci - c^k '^-k) (57) 
[clcU,clcU = (58) 

From the latter it follows that [B^"'\B^^^] = [D^"'\D^^^] = [B^"'\D^^^] = 0, so [A, C] = [Bi, Bj] = 
[Di,Dj]^[Bi,Dj]^0. 

[clck,4cli\ = 2 5k,iclc_k, [cfc4,cfcL,] = -2 4_icLfe Cfc (59) 
[clck,c\ci] = [40^,^4] =0 (60) 

The mode operators B^""^ also satisfy [S^"*) , S^'^^t] ^ -2zi"'\ [Z^"'\ B^"'^^ = B^'^^^ where 

^Mt ^ si™) and 5^"^ = 4""^. See also (31-33). Also, [D^"^), i?('")t] = -2Y:,^"'\ 

[yM 2;(m)t] ^ ^(m)t rpj^^g^ f^j. ^^^.j^ ^j^g ^M^g ^^^^ j^{m),g j^g^ ^ two-body Operator, quasi- 

spin, and single-mode Hamiltonian structure, analogous to those at the A's and C's level. 

Some of the commutators vanish at the A's and C's level because [B^\ S^™-*]!;^^ = 
[D^^lD^i^^ = [B^\Dth = for a,b = 1,2,3. Thus, [A,C^ = [B.^D}] = [A, Dj] = 
[A, A] = = [C,Bi] = 0, and [Za3,Yc3] = [Z^Bi.Yc^] = {Za^^Y^d^ = [Z^Bi.Y^Dr] = 0. 

Also, [Z^3,Ct] = [Z3B^,C1 = [^A3, A^] = [ZzBi,D]\ = 0, [Fcs, ^t] = [y3,3i,At] = [yc3,Sl] = 

[YsDj, i?!] = , and their adjoints vanish. 
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